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ABSTRACT 

We report a possible detection of TeV gamma-rays from the Galactic Center by 
the Whipple 10m gamma-ray telescope. Twenty-six hours of data were taken over an 
extended period from 1995 through 2003 resulting in a total significance of 3.7 standard 
deviations. The measured excess corresponds to an integral flux of 1.6 x 10~ 8 ± 0.5 x 
10~ 8 (stat) ±0.3 x 10~ 8 (sys) photons mT 2 s _1 above an energy of 2.8 TeV, roughly 40% 
of the flux from the Crab Nebula at this energy. The 95% confidence region has an 
angular extent of about 15 arcmin and includes the position of Sgr A*. The detection 
is consistent with a point source and shows no evidence for variability. 

1. Introduction 

The central region of our galaxy is now thought to contain a super-massive black-hole of 
2.6x10 6 Mq (Ghez et al. 2002; Schodel et al. 2002) coincident with the unresolved radio source 
Sgr A* (Balick & Brown 1974). Chandra observations reveal X-ray emission from an unresolved 
point source as well as an extended structure (~ 1.5 arcsec), both of which appear to be physically 
associated with Sgr A* (e.g., Baganoff et al. (2003)). The recent discovery of hour-scale X-ray 
(Baganoff et al. 2001) and rapid IR flaring (Ghez et al. 2004) point to an active nucleus, albeit with 
very low bolometric luminosity compared with the luminosity inferred from the Bondi accretion rate 
or with that which is typical of more powerful AGNs. More recently, INTEGRAL {the International 
Gamma-Ray Astrophysics Laboratory) has detected time-variable 20-100 keV emission from within 
0.'9 of Sgr A*(Belanger et al. 2004). Polarization measurements show the signature of synchrotron 
radiation in a Keplarian accretion disk (Liu & Melia 2002). Taken together, these multi- wavelength 
data are not easily described by a one-component model, and the current theoretical framework 
combines thermal emission from a radiatively inefficient Keplarian accretion flow with synchrotron 
inverse-Compton emission produced either by electrons accelerated in the disk or further out in 
a hypothetical jet-like outflow (e.g. Liu k Melia (2002); Yuan, Markoff, & Falcke (2002); Yuan, 
Quataert, & Narayan (2003)). From the present measurements, the maximum energy of the non- 
thermal electron distribution in the jet models is ambiguous and theories alternately explain the 
high-energy emission as inverse-Compton or the high-energy extension of the synchrotron spectrum; 
gamma-ray measurements may eventually break this degeneracy. 

The EGRET experiment detected a strong unidentified source of GeV gamma-rays marginally 
consistent with the position of the Galactic Center (Hartman et al. 1999). Both the Whipple and 
Cangaroo groups have presented preliminary evidence for TeV emission at the position of Sgr A* 
as well (Buckley et. al. 1997; Tsuchiya et al. 2003; Kosack et. al. 2003). Hooper & Dingus 
(2002) re-analyzed the higher energy gamma-ray data from EGRET and found that the most likely 
position of the EGRET source may be offset from Sgr A*. However, systematic uncertainties in the 
gamma-ray background models and limited angular resolution make the analysis of the source in 
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the Galactic Center region difficult. Observations of the Galactic Center are complicated since Sgr 
A* is surrounded by a dense cluster of stars and stellar remnants (including low-mass X-ray binaries 
and black hole candidates), molecular clouds, and a large structure that may be the remnant of a 
powerful supernova remnant, Sgr A East (Fatuzzo & Melia 2003). Source confusion is particularly 
difficult for high-energy gamma-ray observations given the limited angular resolution of present 
experiments. 

High energy gamma-ray observations of the Galactic Center are also the subject of particular 
theoretical interest given the possibility of detecting halo dark matter in our galaxy (e.g. Bergstrom, 
Ullio, & Buckley (1998)). Sgr A*, at the dynamical center of our galaxy, may well be surrounded 
by a cusp or spike in the dark matter halo distribution (e.g. Dubinski & Carlberg (1991); Navarro, 
Frenk, & White (1996); Gondolo & Silk (1999); Merritt (2003)). Annihilation of these hypothetical 
weakly interacting massive particles could also contribute to the luminosity in the vicinity of Sgr A* 
in the radio through gamma-ray waveband. Annihilation of dark matter would be enhanced by a 
factor proportional to the density squared, and might result in an observable gamma-ray line (from 
direct annihilation to gamma rays) as well as continuum emission (from secondary products of 
annihilation to quarks and fermions) (Silk & Bloemen 1987; Bergstrom 1989; Rudaz 1989; Giudice 
&; Griest 1989; Stecker & Tylka 1989; Jungman & Kamionkowski 1995). The presence of a massive 
black hole could further steepen the density profile of the dark matter halo, producing very high 
radio and gamma-ray fluxes that exceed the observational upper bounds (Gondolo & Silk 1999). 
The details of the halo model on scales < lOOpc and the formation history of the central black hole 
are critical to predicting the gamma-ray flux, but are, unfortunately, still poorly understood. 

Given the limited angular resolution of GeV and TeV instruments, a number of different sources 
could contribute to a signal near the Galactic Center. The key to distinguishing between all of the 
possible emission scenarios is to measure the position, angular extent, variability and spectrum of 
the gamma-ray signal. Here we present first results from an analysis of Whipple telescope data. In 
§2 we describe the observational method and data analysis procedure used to observe the Galactic 
Center at TeV energies. In §3 and §4, we discuss a possible weak detection and consider its impact 
on various gamma-ray production scenarios in §5. 



Season(s) 


F.O.V (°) 


Pixels 


Pixel Diameter (°) 


1995-1996 


3.1 


109 


0.26 


1996-1997 


3.1 


151 


0.26 


1999-2003 


2.4 


379 


0.12 



Table 1: Camera specifications during seasons where Galactic Center data were taken. 
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2. Method 

Imaging Atmospheric Cerenkov Telescopes (IACTs), like the Whipple Gamma-ray Observa- 
tory, detect high energy photons by imaging the flashes of Cerenkov light emitted by secondary 
particles in gamma-ray induced air showers. The Whipple Telescope's 10m mirror focuses the faint 
UV/blue Cerenkov flashes into a camera consisting of 379 photomultiplier tube pixels (for the latest 
seasons — see Table 1 for previous seasons). Off-line software analysis characterizes each candidate 
shower image, separates signal (gamma-ray-like) from background (cosmic-ray-like) events, and 
determines the point of origin and energy of each gamma ray. 

Whipple gamma-ray data are traditionally taken as a series of 28 minute exposures, each of 
which is followed by an off-source run which is offset 30 minutes in right ascension for background 
subtraction. In the case of Sgr A*, data were taken off-source before the on-source observations 
due to a bright star field in the region 30 minutes past the Galactic Center's position. The analysis 
method used here (see §3) was modified to give sensitivity at large zenith angles and to provide a 
two-dimensional map of the TeV emission in a three-degree diameter field-of-view surrounding the 
Galactic Center. 

3. Data Analysis 

The Galactic Center (a = 17 h 45 m 40 s , S = -29°00'28", J2000) transits at a very large zenith 
angle (61°) as seen from the Whipple Observatory (31.7° N latitude) which significantly alters the 
shower geometry and threshold energy. To properly account for the effects of LZA observations, 
special techniques that go beyond the standard Whipple analysis were required. Furthermore, the 
brightness of the galactic plane near the Galactic Center results in an increase in energy threshold 
and, if not compensated for, a systematic bias in the observed excess. Pedestal events, containing 
no image, are injected at random intervals throughout the run for calibration of Poisson fluctuations 
in the night sky background. Both on and off-source data are analyzed in the same manner, and 
Gaussian deviates are added to the pixel signals to bring the background noise up to the same 
level in both runs. After this procedure, only pixels with signals well above the noise level are 
included in further image processing. This Gaussian padding combined with a high software trigger 
threshold largely removes systematic biases arising from brightness differences, but increases the 
energy threshold. 

Using the techniques based on the moment fitting procedure outlined in Reynolds et al. (1993), 
we parameterize the roughly elliptical gamma-ray images by calculating moments of the light distri- 
bution in the camera. Geometric selection criteria based on these parameters allow for the rejection 
of background (e.g. cosmic-ray induced showers). The first moments give the centroid of the im- 
age, the second moments give the WIDTH (minor axis) and LENGTH (major axis) and orientation 
angle of the image. The elongation of the ellipse is used to determine the point of origin of each 
shower image using the formula S = e (l — L ^^3^ ) , where 5 is the displacement of the point of 
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origin from the image centroid and e is the elongation factor (determined by simulations). Of the 
two possible points of origin for each image, the asymmetry (or skew) of the shower is used to 
select the correct one whenever possible. This procedure is similar to that described in Lessard et 
al. (2001). 

To determine the pointing error in the telescope, we look at the pedestal variation for each 
tube. The presence of visible light from a star or other source of sky brightness in the field of 
view adds to the pedestal variance in the corresponding photomultiplier tube. Using this effect, we 
can generate a crude optical image of the sky by de-rotating the camera to a common orientation 
(since the field of view rotates with time for an altitude-azimuth telescope) and accumulating the 
pedestal variations of each pixel into a two-dimensional histogram. Using this technique, an optical 
sky-brightness map is generated for each observation. By comparing the bright spots in the Sgr 
A* image (or special runs where the telescope is pointed at a nearby pointing -check star) we can 
obtain an absolute measure of the pointing error. In addition, we cross-correlate each pair of maps 
to determine the relative pointing offset between them. We only keep runs that have the correct 
star field and have a relative pointing offset (compared with the other Sgr A* runs) which is less 
than the diameter of one pixel. This lessens the possibility of accidentally including an observation 
with large pointing errors or which was mislabeled. Using this method, we find residual pointing 
errors of ±0.1° and an absolute offset of 0.14°. These errors are attributed in part to the fact 
that observations have been made near the balance point of the telescope and near the horizontal 
position where flexure of the optical support structure is maximal. 

Sgr A* transits at roughly 30° elevation as observed from the latitude of the Whipple Observa- 
tory. Large zenith angle (LZA) observations require several modifications to the standard analysis 
due to changes in shower geometry. In addition, we have data on Sgr A* that spans several epochs 
during which the Whipple camera was upgraded twice (see Table 1). In order to combine all of 
these data, we needed our gamma-ray selection criteria to scale properly with changes in zenith 
angle, the camera throughput factor (i.e. the ratio of photons hitting the mirror to digital counts 
in the digitized PMT signals), and the pixel size. In addition, for subsequent spectral analysis, we 
designed our gamma-ray selection procedure to scale with energy to minimize spectral biases. 

To achieve this goal, we started with the standard Whipple Telescope data analysis procedure 
(Lessard et al. 2001) and developed a new set of gamma-ray selection criteria that scale with zenith 
angle and energy according to a semi-empirical form derived from simulations and optimization of 
LZA Crab Nebula data. Our criteria also incorporate the geometric effects of differing cameras, 
allowing us easily to combine all of the data present for the Galactic Center. The Crab Nebula 
was used for optimization and calibration because it is a bright, steady gamma-ray source with a 
known spectrum. In addition, an independent set of Crab Nebula observations taken at LZA were 
used to verify that the selection criteria were more efficient than the standard analysis procedure. 

To correct for changes in the overall light sensitivity (throughput) of the camera between 
epochs, we look at muon events in data taken at each epoch. Muons, which show up as bright arcs 
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of Cerenkov light in the camera, are useful for throughput calibration because the light per unit 
arc-length from a muon event is nearly constant regardless of the impact parameter and angle of 
the trajectory. A measure of the throughput of the telescope can be found by making a histogram 
of the signal/ arc- length distribution of a large set of muon events. To get an absolute calibration, 
this distribution for real observations was compared with a set of simulated muon events which 
were generated using the Grinnell/ISU (GrISU) simulation package to produce a large number of 
H and He showers that in turn produce muons as secondary particles. For instance, the relative 
throughput was found to have changed by a factor of 2.22 between 1995 and 2001. We scale the 
software trigger threshold and energy estimator by this factor so the trigger cuts are consistent 
across observations. This method also serves to calibrate the simulations used to determine the 
peak energy at LZA (see §4). 

As one observes at increasing zenith angles, the distance to the core of the air-shower increases 
and thus the angular size of the shower and parallactic displacement of the image centroid are 
reduced. To derive the scaling laws, we first assume that the WIDTH and LENGTH of gamma-ray 
air shower images is approximately proportional to cos a 0, where 9 is the zenith angle and a is a 
constant. Additionally, air-shower simulations show that LENGTH and WIDTH also scale as the 
logarithm of the energy, which is proportional to the SIZE (total camera signal) of the event. 

Combining these results, and removing the effects of the finite pixel size of the camera (cr p i x ) 
and the point spread function of the telescope (cr ps f), the measured LENGTH (L), WIDTH (W), 
and the distance to the image centroid (D) can be converted to scaled values L', W' and D' by the 
following equations: 



E = ln(SIZE) - 8.0 



W' ~ 



T 2 _ 2 -2 
^ a pix _ a psf 

cos 1 - 5 9 
W 2 -a 2 . -a 2 r 

w pix w psf 

cos 1 - 2 9 



D' 



0.023S 

0.020S 
D 



cos 9 



(1) 

(2) 
(3) 



The constant factors and cosine powers were derived from simulation fits and by optimization on 
Crab Nebula data (taken at a range of zenith angles) and then compared with simulations for 
calibration. Data selection criteria 0.125° < V < 0.3°, 0.05° <W'< 0.135°, 0.28° < D' < 2.2° 
are applied to select candidate gamma-ray events (examples of these cuts are shown graphically in 
Figure 1). Cuts based on these intrinsic parameters were verified to be independent of zenith angle 
and camera design by application of this method to independent Crab Nebula data taken over the 
period 1994-2003. 
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Fig. 1. — Plots of our selection criteria for gamma-ray-like events for two of the parameters 
(LENGTH and WIDTH) for two different zenith angles (0° and 60°). The dotted lines show 
the upper and lower cuts on the respective parameters as a function of the SIZE (total signal) 
parameter. The dots are from simulations of gamma rays with a range of energies. 



4. Results 



We have combined observations of Sgr A* from 1995 through 2003 resulting in 26 hours of on- 
source exposure at an average zenith angle of 61°. To determine the pointing offset, observations 




Fig. 2. — Optical image of a star (Sgr 72) used to check the telescope's pointing at low elevation. 
This image shows that the telescope has an offset of 0.14 degrees down and to the right of camera 
center position (0,0). 
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Fig. 3. — A gamma-ray image of the region around Sgr A*. The image is of excess counts with 
overlaid significance contours (1 standard deviation per contour). The axes are labeled in degrees 
from the assumed camera center. The true center position of the camera, which is not exactly at 
(0,0) due to flexing of the telescope at low elevation, is marked with a cross. The dashed lines are 
the RA and Dec contours at this position. Also shown (as a light contour) is the 99% confidence 
region for the EGRET observations (Hooper & Dingus 2002). 

were taken centered on a nearby bright star (Sgr 72) which is at the same elevation as Sgr A*. 
Using the sky brightness map technique outlined earlier, we determined that the telescope had a 
pointing offset of 0.14°, as shown in Figure 2. Figure 3 shows the resulting 2-D map of gamma-ray 
excess with overlaid significance contours. The true center of the camera, correcting for the offset, 
is plotted as a cross in the image. This image shows a 4.2 standard deviation (<r) excess at the 
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Fig. 4. — A gamma-ray image of the Crab Nebula taken at large zenith angle (~ 62°) using the 
same analysis procedure used for Sgr A*. The offset and pointing variations can be seen in the 
resulting image. 



corrected center position. To check the robustness of this result, we re-ran the analysis ten times 
to account for variations due to the Gaussian padding. We find the average significance at the 
corrected center position is (3.7 ± 0.13) a, somewhat below the initial result. For reference, in 
Figure 4 we show the results of the same analysis procedure applied to 10 hours of observations 
of the Crab Nebula at a similar zenith angle range. Note that the significance of 6. la of the Crab 
detection at the offset position is substantially higher than the result of 3.8a obtained applying the 
standard small zenith angle analysis procedure to these LZA data. Also, the similar angular extent 
in the two results indicates consistency with a point source within a 99% confidence region of radius 
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ps 15 arcmin. A more detailed analysis of the angular extent would require a prior hypothesis about 
the nature of the source. 

To determine the peak energy of the detected flux from Sgr A*, we simulated gamma rays with 
a Crab Nebula spectrum (with integral spectral index 7 = 1.58) (Mohanty et al. 1998; Aharonian et 
al. 2000) and a zenith angle of 61°, and analyzed the resulting data with a detector simulation and 
our analysis software. We determined the peak detected energy to be 2.8 TeV. Noise padding 
and a higher trigger threshold makes this result slightly higher than has been previously reported 
(Krennrich et al. 1999), and we estimate a 20% systematic error in this energy threshold. We then 
analyzed a set of real LZA Crab Nebula data runs to find the Crab count rate and compared this 
to the corresponding rate for Sgr A*. The integral flux for Sgr A*, normalized to the Crab, is then: 

F SgrA (> 2.8 TeV) = A , Crab • &^XL! . ^ (4) 

7 ^Crab 

Where N^crab is the flux normalization factor for the Crab Nebula (3.12 x 10 _7 m _2 s _1 ), 7 is the 
integral Crab spectral index, and i?sgrA* and -Rcrab are the corresponding Sgr A* and Crab Nebula 
gamma-ray count rates. From the LZA Crab data, we find a gamma-ray rate of i?crab(> 2.8 TeV) = 
0.501 ± 0.087 photons min -1 and from Sgr A* we obtain an average rate of i?sgrA*(> 2.8 TeV) = 
0.205 ± 0.057 photons min -1 . Hence, the gamma-ray flux from the Galactic Center region above 
2.8 TeV is 1.6 x 10" 8 ± 0.5 x 10~ 8 (stat) ± 0.3 x 10" 8 (sys) photons mT 2 s _1 , or about 0.4 times that 
of the Crab Nebula (the flux error includes the uncertainty in the Crab Nebula measurement). 

The count rate from Sgr A* is shown as a function of time in Figure 5. To determine the 
probability for steady emission, a x 2 fit of a constant function (f(t) = A) was applied to this data 
and, for comparison, to a series of data taken of Mrk 421 (a source which is known to be highly 
variable) at a similar zenith angle range as Sgr A*. The total significance of this Mrk421 data 
sample was 2.3a. The Sgr A* data yields a constant count rate of 6.12 ± 1.59 7 • min -1 with a 
reduced \ 2 of 1.13 (with 54 degrees of freedom), which corresponds to a 25% probability that there 
is no variability. The result for Mrk 421 yields a constant count rate of 6.86 ± 6.13 7 • min -1 with 
a reduced \ 2 of 3.03 (with 6 degrees of freedom) and a 1.2% chance of no variability. 



5. Discussion 

The TeV excess observed near the position of the Galactic Center is unlikely to have occurred 
by chance and constitutes a probable, as yet unconfirmed, detection of a new TeV source. Possible 
systematics that could contribute to a false detection include the effects of additional noise from the 
relatively bright off-source region. While we have largely corrected for these effects, some systematic 
uncertainties remain. We have taken into account trials factors by formulating an explicit a priori 
hypothesis that we would only look for emission at the exact position of the Galactic Center after a 
pointing correction was applied. Statistical variations in the analysis method (due to the addition of 
simulated noise in padding) have been taken into account by repeating the analysis ten times, and 
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Fig. 5. — Flux of Sgr A* as a function of time. Each data point represents a single 28 minute run. 
Time gaps in the data have been removed where indicated. The dashed line is a least squares fit 
of a constant function to the data. 

taking the average significance, giving a conservative estimate of 3.7a for the detection significance. 

The lack of significant variability in our data makes it difficult to uniquely identify the source 
with a compact point source such as Sgr A*, but inspires some confidence in the stability of our 
observations at large zenith angle. Note that the analysis procedure was designed to mitigate 
against changes in the count rate due to variations in the instrument. The same ISU simulation 
package was used here to analyze Whipple observations of the Crab Nebula giving a spectrum in 
good agreement which that measured a small zenith angle (Krennrich et al. 1999). In the past, our 
group reported a positive excess of 2. 4a for 1995-1997 observations (Buckley et. al. 1997) and 
2.4a for 1999-2003 observations (Kosack et. al. 2003) at the position of Sgr A*. The combined 
significance for our refined results is consistent with these earlier analyses. The large error circles 
for both EGRET (7.2 arcmin) and Whipple (15 arcmin) observations make identification with a 
particular source difficult, but given the dearth of TeV sources, an accidental angular coincidence 
of a new source along the line of sight is unlikely, and it is probable that the emission comes from 
a non-thermal source physically near the Galactic Center. 

The high level of emission « 0.4 Crab at a distance of roughly four times that of the Crab 
Nebula, qualifies this as an unusually luminous galactic source. Previous TeV observations of 
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relatively nearby galactic sources such as X-ray binaries and shell-type and plerionic supernovae 
have produced numerous upper-limits, or (at best) unconfirmed detections, making the detection 
of such an object at 8.5kpc even more unlikely. Mayer-Hasselwander et al. (1998) came to a similar 
conclusion about the GeV emission based on the high luminosity of the EGRET unidentified source 
and lack of significant variability. If the Sgr A East supernova shock were the source of the EGRET 
gamma rays, it would have been an unusually intense explosion (Khokhlov & Melia 1996) and a 
density of 1000 cm -3 and magnetic field of B ~ 0.18mG (well above the canonical values) would be 
required (Fatuzzo & Melia 2003). While a typical galactic source such as an SNR, pulsar, or stellar 
mass black hole is unlikely, an active nucleus at our Galactic Center is still a viable possibility and 
the detection of correlated variability in future gamma-ray and X-ray observations is required to 
make the identification with Sgr A* compelling. If we associate this emission with either the super- 
massive black hole Sgr A* or supernova remnant Sgr A East, the observed emission could come 
from self-Compton scattering by electrons with energies up to at least 2.8 TeV or from pion-decay 
gamma-rays from primary protons of even higher energy (kinematics require their energy to be at 
least several times the maximum gamma-ray energy). 

The lack of significant variability and the consistency with the Galactic Center position allow 
more exotic possibilities such as the annihilation of very high mass (> 2 TeV) dark matter particles 
at the Galactic Center. While not particularly constraining, these results and a more detailed 
analysis of the spectrum will be used to derive upper limits for dark matter annihilation in a 
subsequent paper. 

In summary, based on 26 hours of data taken with the Whipple 10m gamma-ray observatory, we 
report a probable detection of a gamma-ray point source consistent with the position of the Galactic 
Center. We also describe a modified analysis procedure that we have developed to analyze gamma- 
ray data taken by Imaging Atmospheric Cerenkov telescopes operating at large zenith angles. The 
continued observation of TeV gamma rays from the Galactic Center has important theoretical 
implications for understanding a variety of astrophysical phenomena in the nuclear region of the 
Milky Way. 
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